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Abstract 

The formulas which describe the time evolution of radiation produced by a rela- 
tivistic particle moving in a crystal are derived. It is shown that the conditions are 
realizable under which parametric (quasi-Cherenkov) radiation, transition radiation, 
diffracted radiation of the oscillator, surface quasi-Cherenkov and Smith-Purcell ra- 
diation last considerably longer than the time r p of the particle flight through the 
crystal. The results of carried out experiments demonstrate the presence of addi- 
tional radiation peak appearing after the electron beam has left the photonic crystal. 



Introduction 



At present, the processes of diffracted radiation of photons by relativistic parti- 
cles passing through crystals (natural or artificial spatially periodic structures) 
are intensively studied both theoretically and experimentally. Worthy of men- 
tion are such types of diffracted radiation as parametric (quasi-Cherenkov) 
radiation and diffracted radiation of a relativistic oscillator [IH3] . It should be 
noted, however, that until now, theoretical and experimental analysis of radi- 
ation produced by a relativistic particle passing through a crystal has focused 
on spectral-angular characteristics of radiation. Nevertheless, it was shown 
in [5l|6] that because of diffraction, photons produced through radiation in 
crystals have group velocity v? r , which is appreciably smaller than the veloc- 
ity v of a relativistic particle. As a result, the situation is possible in which 
radiation from the crystal still continues after the particle has passed through 
it [S1E]- This enables studying time evolution of the process of photon radi- 
ation produced during the particle transmission through the crystal (natural 
or photonic), or during the particle flight along the surface of such crystals. In 
the present paper the formulas are derived, which describe the time evolution 
of radiation produced by a relativistic particle moving in a crystal. It is shown 
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that the conditions are realizable under which parametric (quasi-Cherenkov) 
radiation, transition radiation, diffracted radiation of the oscillator, surface 
quasi-Cherenkov and Smith-Purcell radiation last considerably longer than 
the time r p of the particle flight through the crystal, i.e., much longer than 

T p < 10" 9 S. 



1 Spectral-angular distribution of radiation produced by a particle 
transmitted through a crystal 



Let us first recall the conventional consideration of the radiation process in 
crystals [HE]. 

Both the spectral-angular density of radiation energy per unit solid angle Wa u 
and the differential number of emitted photons dN^ui = 1/ftio ■ can be 
easily obtained if the field E(f,co) produced by a particle at a large distance 
r from the crystal is known [3] 



er 
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E(f,u) 



The vinculum here means averaging over all possible states of the radiating 
system. In order to obtain E(r,u), Maxwell's equation describing the interac- 
tion of particles with the medium should be solved. The transverse solution 
can be found with the help of Green's function of this equation, which satisfies 
the expression: 

G = G + G ^(e-1)G, (2) 



Go is the transverse Green's function of Maxwell's equation at e — 1. It is 
given, for example, in [19] . 

Using G, we can find the field we are concerned with 

E n (f,u) = / G n i(ry,u) l —j i(r,ou)d?r', (3) 
J c 



where n,l — x, y, z, joi{r, u) is the Fourier transformation of the e-th compo- 
nent of the current produced by a moving beam of charged particles (in the 
linear field approximation, the current is determined by the velocity and the 
trajectory of a particle, which are obtained from the equation of particle mo- 
tion in the external field, by neglecting the influence of the radiation field on 
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the particle motion). Under the quantum-mechanical consideration the cur- 
rent jo should be considered as the current of transition of the particle-medium 
system from one state to another. 

According to [21 E], Green's function is expressed at r — > oo through the 
solution of homogeneous Maxwell's equations E^~'(r,ui) containing incoming 
spherical waves: 



hmGUr>->) = — E<47 ( r ~>)' ( 4 ) 
r s 

r — > oo 

where (? is the unit polarization vector, s — 1, 2, e 1 _L e 2 _L k. 

If the electromagnetic wave is incident on a crystal of finite size, then at r — > oo 



ikr 

E { ~ )s {r,u) = eV fer ~+ const 



and one can show that the relation between the solution E k and the solution 
of Maxwell's equation E^ +S> {k^u) describing scattering of a plane wave by the 
target (crystal), is given by: 

E<f>* = E {+ J s (5) 

k — k v ' 



Using ©, we obtain 



e ikr iuj 



r c 2 s 



E n (r,u) = — / E { t )s \f^)UrMd"r'. (6) 



As a result, the spectral energy density of photons with polarization s can be 
written in the form: 
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to 



An 2 c 2 



E 



(-)s* 



r, u)jo(r, u)d 3 r 



(7) 



j (r,a;) = / e wt j (r,u)dt = eQ / e^v(t)5(f - f(t))dt 



where eQ is the charge of the particle, v(t) and f(t) are the velocity and the 
trajectory of the particle at moment t. By introducing ([8]) into (J7j) we get 
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(9) 
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Integration in (Q is carried out over the whole interval of the particle motion. 
It should be noted that the application of the solution of a homogeneous 
Maxwell's equation instead of the inhomogeneous one essentially simplifies the 
analysis of the radiation problem and enables one to consider various cases of 
radiation emission taking into account multiple scattering. 

Using equations ©-Q, one can easily obtain the explicit expression for the 
radiation intensity and that for the effect of multiple scattering on the process 
under study |3j[8j|9]. 

Consider, for example, the PXR radiation. Let a particle moving with a uni- 
form velocity be incident on a crystal plate with the thickness L being L <C L c , 
where L c = {uq)~ 1 ^ 2 is the coherent length of bremsstrahlung q = 9 /4 and 9 
is the mean square angle of multiple scattering. The latter requirement allows 
neglecting the multiple scattering of particles by atoms. A theoretical method 
describing multiple scattering effect on the radiation process is given in [ID] . 



According to ([9]), in order to determine the number of quanta emitted by a 
particle passing through the crystal plate, one should first find the explicit 
expressions for the solutions E^ s . As was mentioned above, the field Ejf 

can be found from the relation E^ s = (E^ s )* if one knows the solution 

E^ s describing the photon scattering by the crystal. 

In the case of two strong waves excited under diffraction (the so-called two- 
beam diffraction case [TT]). one can obtain the following set of equations for 
determining the wave amplitudes (see [12]): 



(£ - 1 - a) £r«;4~" = o- <«» 

Here k? = k + f, f is the reciprocal lattice vector, xo-> Xr are the Fourier 
components of the crystal susceptibility. It is well known that the crystal is 
described by a periodic susceptibility (see, for example, [11] : 

Xif) = J2 Xr exp(ifr). (11) 



c s = e^e^, where (^(ef) are the unit polarization vectors of the incident and 
diffracted waves, respectively. 

The condition for the linear system (11 01) to be solvable leads to a dispersion 
equation that determines the possible wave vectors A; in a crystal. These wave 
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vectors are convenient to present in the form: 



k 



/IS 



k + k; s n, 



where /i = 1, 2; A" is the unit vector of a normal to the entrance crystal surface 
which is directed into the crystal, 



ei(2) s 



t [(1 + /3i)Xo - ± 7 {[(1 - ft)*) + Pia B f 



+A(3 1 C 2 sX rX-^} 



-1/2 



(12) 



«e = (2A;r + r 2 )A; 2 is the off-Bragg parameter («g = if the exact Bragg 
condition of diffraction is fulfilled), 



70 = n y ■ N, n 7 = -, p x = — , 71 = n 7T • AT, n 7T 
k 7i 



k + T 

\k + f\ 



The general solution of (TlOT) inside a crystal is: 
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E [eM M exp(i& MS r) + e^A TA1 exp(zA; A1ST f) 



(13) 



Associating these solutions with the solutions of Maxwell's equations for the 
vacuum area, one can find the explicit form of Ei (r) throughout the space. 
It is possible to discriminate several types of diffraction geometries, namely, 
the Laue (a) and the Bragg (b) schemes are most well known. 

(a) Let us consider the PXR in the Laue case. 



In this case, the electromagnetic waves emitted by a particle in both the 
forward and the diffracted directions leave the crystal through the same surface 
(k z > 0, k z + t z > 0), the z-axis is parallel to the normal A" (where A" is the 
normal to the crystal surface being directed inside a crystal). By matching the 
solutions of Maxwell's equations on the crystal surfaces with the help of ffTU]) . 
f[T2l . (|13p . one can obtain the following expressions for the Laue case: 
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x9(L-z)9(z) + S s e ik %z-L), (14) 
where 



S1.2s 




9(z) 



Substitution of (fT4"|) into © gives for the Laue case the differential number of 
quanta of the forward directed parametric X-rays with the polarization vector 
e s : 



e 2 Q 2 w 2 
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(15) 



where T = L/cjq is the particle time of flight; e\ || [fcr]; e 2 || [kei] 



One can see that formula ( JTBT) looks like the formula which describes the 
spectral and angular distribution of the Cherenkov and transition radiations 
in the matter with the index of refraction = k Zfls /k z = 1 + n^ s /k z . 

The spectral angular distribution for photons in the diffraction direction k T = 
k + t can be obtained from (TT5"|) by a simple substitution 



C s t C ST , 
£1(2)8 = =t 

k y Jv-r , 



XtC 



/LtS ' 
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k/j, s + r - 



(b) Now let us consider PXR in the Bragg case. In this case, side by side with 
the electromagnetic wave emitted in the forward direction, the electromagnetic 
wave emitted by a charged particle in the diffracted direction and leaving 
the crystal through the surface of the particle entrance can be observed. By 
matching the solutions of Maxwell's equations on the crystal surface with the 
help of ( TTU1) . (TT2l) . (fl3l) . one can get the formulas for the Bragg diffraction 
schemes. 
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It is interesting that the spectral angular distribution for photons emitted in 
the forward direction can be obtained from (|T5|) by the following substitution, 

7i%. = ^ )s -Xo (16) 



The spectral angular distribution of photons emitted in the diffracted direction 
can be obtained from f[T5"j) by substitution 



where 



T PlCsXr (17\ 

7l(2)s , , , v J£( E 2(1) s _£ l(2) s) L ' ^ ' 



Let us note that the above formulas fully describe parametric ( quasi- Cherenkov) 
radiation in natural and photonic crystals and they certainly include that con- 
tribution to radiation, which goes over to ordinary transition radiation, if the 
radiation is considered outside the region of diffraction reflection. A descrip- 
tion of diffracted radiation of a relativistic oscillator is given in [TJ[2] and the 
reference therein. 

Let us take notice of the fact that in photonic crystals built from metal threads 
with the diameter smaller than or comparable with A, the value of x( r ) is 
practically independent on r. As a result, it is possible to effectively excite 
radiation in, e.g., the terahetrz range in a lattice with a period of several 
millimeters. 

When a particle travels in a vacuum near the surface of a spatially periodic 
medium, new kinds of radiation arise [131 E] ~ surface parametric (quasi- 
Cherenkov) X-ray radiation (SPXR) and surface DRO (see Figure [[)). This 
phenomenon takes place under the condition of uncoplanar surface diffraction, 
first considered in [15] . 

The solution of Maxwell's equation (r ) in this case of uncoplanar surface 
diffraction was obtained in [15]. It was shown that the surface diffraction in 
the two-wave case is characterized by two angles of total reflection (several 
angles in the case of multi-wave diffraction [IS]). The solution obtained in JTH] 
contains the component, which describes the state that damps with growing 
distance from the surface of the medium, both within the material and in the 
vacuum, and which describes a surface wave, i.e., a wave in which the energy 
flux is directed along the boundary of the surface of a spatially periodic target 
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photonic crystal 

Figure 1. Surface diffraction of a radiated photon 

(see review [17]). According to [15], this solution, which describes scattering 
of a plane wave by the target under the surface diffraction geometry, can be 
written in the form: 



E 



e s e ikf + A 8 (k, u)e lklf + B s (k, u)e 



ik\r 



ik2r 



;is) 



where the wave vector in a vacuum k = (kt,k±), k\ = (kt,—k±), k 2 = 
(k 2 t, —k2±)i \k 2 ±\ = \Jk 2 — k% t , k 2t — k t + r, k t is the component of the wave 
vector that is parallel to the surface, f is the reciprocal lattice vector, u is 
the photon frequency The amplitudes A s and B s are given in [SlIH]- Substi- 

distribution of SPXR and DRO^ 



tuting the solution E 



into (j3J), we can find the spectral-angular 



2 Time dependence of the intensity of radiation produced by a 
particle transmitted through a crystal 

The intensity I(t) of radiation produced by a particle which has passed through 
a crystal can be found with known intensity of the electric field E(r, t)) (mag- 
netic field H(f,t)) of the electromagnetic wave, which is produced by this 
particle [T8], 

lit) = — \E(rM 2 r 2 dn, (19) 
Atx 

where r is the distance from the crystal, which is assumed to be larger than 
the crystal size. 

The field E(f,t) can be presented as an expansion in a Fourier series 

E(f,t) = — [ E(f,u)e- iuJt du. (20) 

Z7T J 
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According to the results obtained in [3|TT|l9] , at a long distance from the crystal, 
the Fourier component can be written as follows: 



E{f,t) 



e lkr iui 



r c* 



s 



E 



V. 



(21) 



where i — 1, 2, 3 (and correspond(s) to the coordinate axes x, y, z), e\ is the i- 
component of the wave polarization vector e 8 ; s = 1, 2; E~~^ a is the solution of 
Maxwell's equations describing scattering of a plane wave with a wave vector 
k = and the asymptotic of a converging spherical wave, 

j(r» = Jj(f,t)e^dt (22) 



j(f,u) = Qv(t)5(f — r(i)) is the is the current density of the particle with 
charge Q, r{t) is the particle coordinate at time t. 

The explicit form of the expressions E^~' s describing diffraction of the elec- 
tromagnetic wave in a crystal in the Laue and Bragg cases is given in [3j,[8|[T2] 
(See Section 1). 

Now let us take a closer look at the expression for the amplitude A(u>) of the 
emitted wave: 

Af (u;) = ^J (r*, uW, u)d 3 r>. (23) 



Using ( |22i) . ( 1231 can be recast as follows 
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(24) 



Recall that )s * = E ( t )s , where the field E {+ J s is the solution of Maxwell's 

k —k —k 

equations describing scattering by a crystal of a plane wave with wave vector 
(—A;) and the asymptotics of a diverging wave at infinity. According to (I24p . 
the radiation amplitude is determined by the field E^ s taken at point r(t) of 
particle location at time t and integrated over the time of particle motion. 

Let us consider in more detail the constant motion of a particle in passing 
through the crystal. In this case, parametric quasi- Cherenkov radiation can 
appear [HE], which includes, as a particular case, diffracted transition radia- 
tion. The explicit formulas for the radiation amplitude in the case of two- wave 
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diffraction of photons in crystals for the Laue and Bragg geometries are given 
in 1310112] (Section 1). 

From ( 120 p . ( 121 p . and ( 12 3 p follows that the expression for the electromagnetic 
wave emitted by the particle passing through the crystal (natural or photonic) 
can be presented in a form: 

Ei(f, t) = ^- r Y.<j A^e-^-Vdu, (25) 



i.e.,E l (r,t) = lEe s i Al(t-l). 

From ( 1251) follows that the time dependence of the form of the pulse I(f, t) (E(r, t)) 
of radiation generated by a particle passing through the crystal is determined 
by the dependence of the radiation amplitude A s ^(u) on frequency. Accord- 
ing to the explicit expression for the radiation amplitudes given in [3| [8|, [12] . 
the radiation amplitudes A s ^(oj) can be presented as sums proportional to the 
amplitudes of diffraction reflection from the crystal and to the amplitude of 
wave transmission through the crystal. For example, for the case of forward 
parametric radiation in the Laue geometry 



H=l,2 

1 1 r •/ ,r, *- r s^ r, 1 

(26) 



oj — kv oj — (k + K/j, s N)v_ 



Thus, the time dependence of the from of the radiation pulse is determined 
by the time dependence of the radiation amplitude A~(t — ^). 

By way of example, let us consider the characteristics of the time dependence 
of radiation produced by a particle passing through the crystal for a wave 
packet passing through the crystal [2J[5j[6] 

Let us consider the pulse of electromagnetic radiation passing through the 
medium with the index of refraction n(oj). The group velocity of the wave 
packet is as follows: 

fau^y l = c 



where c is the speed of light, oj is the quantum frequency. 
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In the X-ray range ( ~tens of keV) the index of refraction has the universal 
form n(u) = 1 — ^% > is the Langmuir frequency. Additionally, n — 1 ~ 

10~ 6 <C 1. Substituting n(u) into ff27|) . one can obtain that v gr ~ c ^1 — 

It is clear that the group velocity is close to the speed of light. Therefore the 
time delay of the wave packet in a medium is much shorter than the time 
needed for passing the path equal to the target thickness in a vacuum. 

AT = ~ — I < -■ (28) 

v ar c cu z c 



To consider the pulse diffraction in a crystal, one should solve Maxwell's equa- 
tions that describe a pulse passing through a crystal. Maxwell's equations are 
linear, therefore it is convenient to use the Fourier transform in time and to 
rewrite these equations as functions of frequency: 




- i 



where Xijif^) is the spatially periodic tensor of susceptibility; i,j = 1,2,3 
repeated indices imply summation. 

Making the Fourier transformation of these equations in coordinate variables, 
one can derive a set of equations associating the incident and diffracted waves. 
When two strong waves are excited under diffraction (the so-called two-beam 
diffraction case), the following set of equations for determining the wave am- 
plitudes can be obtained: 

(S - 1 - Xo) El - c sX ^Ei r = 

(30) 

(5 - 1 - Xo) El - c sX ,Ei = 



Here k is the wave vector of the incident wave, k? = k + r, f is the reciprocal 
lattice vector; Xo,Xr are the Fourier components of the crystal susceptibility: 

X(r) = Xr exp(zff) (31) 



C s = e s e j f, e s (e/) are the unit polarization vectors of the incident and 
diffracted waves, respectively. 

The solvability condition for the linear system (130]) leads to a dispersion equa- 
tion that determines the possible wave vectors k in a crystal. It is convenient 
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to present these wave vectors as: 



where \i — 1, 2; N is the unit vector of a normal to the entrance surface of the 
crystal, which is directed into the crystal, 

4 1,2) = t[(1 + P)Xo - P<* B ] ± - {[(1 + /3)xo - - 2 Xo f + ^XfX-r} 1 



olb = (2kr + r 2 )A; 2 is the off-Bragg parameter (ag = when the Brag 
condition of diffraction is exactly fulfilled), 

- at - k a ^° jvr - 

7o = n 7 • N, n 7 = -, = — , 71 = n 7T • AT, n 7r - 



k 7i |fc + r| 

The general solution of equations ( 1291) . ( |30j) inside a crystal is: 
2 

if) = X) [ g exp(i^ s r) + e T s A r ^ exp(i£ MST r)] (33) 
(1=1 



Associating these solutions with the solutions of Maxwell's equation for the 
vacuum area one can find the explicit expression for E~(f) throughout the 
space. It is possible to discriminate several types of diffraction geometries, 
namely, the Laue and the Bragg schemes, which are most well-known 



In the case of two-wave dynamical diffraction, the crystal can be described by 
two effective indices of refraction 

W = 1 | Xo (l + f3)-f3a± ^(l-^+HH^M-r} • (34) 



The diffraction is significant in the narrow range near the Bragg frequency, 
therefore Xo an d Xt can be considered as constants and the dependence on 

oj should be taken into account for a = 27r ^( 27r ^+ 2 b. ) = — C 2 ^) (u — 00 b), 
where k = -; 2ix^f is the reciprocal lattice vector which characterizes the set 
of planes where the diffraction occurs; Bragg frequency is determined by the 
condition a = 0. 
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From ( p7|) . (|34|) one can obtain 



„(i.2)« = T 



Ti a/(xo(1-/3)+/3«) 2 +4/3C s XtX- 



In the general case (xo(l _ 0) + /3a) — 2v^X0) therefore the term that is 
added to n^-' 2 \co) in the denominator f )35l) is of the order of 1. Moreover, 
v gr significantly differs from c for the antisymmetric diffraction ^> 1). It 
should be noted that because of the complicated character of the wave field 
in a crystal, one of v g l ) s can appear to be much higher than c and negative. 
When /3 is negative the radicand in f l33|) can become zero (Bragg reflection 
threshold) and v gr — > . It should be noted that in the presence of a variable 
external field, a crystal can be described by the effective indices of refraction 
which depend on the external field frequency Q . Therefore in this case v gr 
appears to be the function of Q . This can be easily observed in the conditions 
of X-ray-acoustic resonance. The performed analysis allows one to conclude 
that the center of the X-ray pulse in a crystal can undergo a significant delay 
AT ^> - available for experimental investigation. Thus, when (3 = 10 3 , I = 0.1 
cm and l/c ~ 3 • 10~ 12 , the delay time can be estimated as AT ~ 3 ■ 10 _9 sec. 

Let us study now the time dependence of the delay law of radiation after 
passing through a crystal. Assuming that B(u) is the reflection or transmission 
amplitude coefficients of a crystal, one can obtain the following expression for 
the pulse form 

E(t) = — / B(u)E (u)e' iult du = I Bit- t')E (t')dt' . (36) 

where Eq(lij) is the amplitude of the electromagnetic wave incident on a crys- 
tal 



In accordance with the general theory, for the Bragg geometry, the amplitude 
of the diffraction-reflected wave for the crystal width much greater than the 
absorbtion length can be written as [22] : 

B s {u) = (37) 
1 



2x T 



Xo(l + v / (Xo(1-|/3|)-|/3|«) 2 -4|/3|C sXt X- 



In the absence of resonance scattering, the parameters Xo an d X±r can be 
considered as constants and frequency dependence is defined by the term 
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ot = — ( - 2 , 7r 3 r ' ) (to — lu b ). So, B s {t) can be found from 



B s (t) = — (38) 

x / {%o(l + \P\ a - v / (xo(l-|/3|)-|/3|«) 2 -4|/3|C sX rX-.} e^du. 



The Fourier transform of the first term results in 8(t) and we can neglect it 
because the delay is described by the second term. The second term can be 
calculated by the methods of the theory of function of complex argument: 

Bs (t) = - — W 2jl ( a ^) e ^+^) { g (i)] (39) 
AXr k%u B t 



or 



B s (t) 



iy/\P\ Ji{a s t) 



-i(oj B +Auj B )t 



ad 



9(t), 



(40) 



where 



j^jW ' UB ~ " |/3|(2vrr)2 ' 



Since Xo an d Xr axe complex, both a s and Au B have real and imaginary 
parts. According to (j3"9|) - (j4"0"{) . in the case of Bragg reflection of a short pulse 
(the pulse frequency bandwidth 3> frequency bandwidth of the total reflection 
range) both the instantly reflected pulse and the pulse with amplitude under- 
going damped beatings appear. Beatings period increases with \/3\ grows and 
Xt decrease. Pulse intensity can be written as 

2 

e- 2lmAuJBt 6(t). (41) 



Lit) ~ \B s (t)\ 2 = l -f 



Ji(a s t) 



at 



It is evident that the reflected pulse intensity depends on the orientation 
of photon polarization vector e s and undergoes the damping oscillations on 
time. 

Let us evaluate the effect. Characteristic values are ImAus ~ ^Xo^b and 
Ima ~ Im *iyB . For 10 keV for the crystal of Si Imxo = 1, 6 ■ 1(T 7 , for LiH 
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ImXo = 7,6 • 10- u ,lmxr = 7 • 10 -11 , for LiF Im^o ~ 10~ 8 . Consequently, 
the characteristic time r for the exponent decay in f HTj) can be estimated as 
follows (oj b = 10 19 ): 

for Si the characteristic time r ~ 10~ 12 sec, for LiF the characteristic time 
t ~ 10~ 10 sec, for LiH the characteristic time r ~ 10 -9 sec!! 

The reflected pulse also undergoes oscillations, the period of which increases 
with growing \/3\ and decreasing Rex T . This period can be estimated for = 
10 2 and Rexr ~ 10~ 6 as T ~ 10" 12 sec (for Si, LiH, LiF). 

When the resolving time of the detecting equipment is greater than the os- 
cillation period, the expression (HTj) should be averaged over the period of 
oscillations. Then, for the time intervals when Kea s t 3> 1, ImAc^t <C 1 the 
delay law (|4T|) has the power function form: 

i.(t) ~ r 3 . 

In the case of multi-wave diffraction, the time delay for the photon exit from 
the crystal will be even more appreciable. 

For an artificial spatially periodic medium (diffraction grating, photonic crys- 
tal), the parameter goi can vary over a wide range. For example, according 
to [23], for a photonic crystal built from tungsten threads of 100/im in diame- 
ter, the parameter g i ~ \ has the value of poi ~ 10 -2 in a 10 GHz range. As 

a result, in this range we have T (10 GHz)~ ~ ■ 10 -9 . At the same 

time, in the terahertz range, T (1 THz) due to the drop of poi (T increases 
proportionally to u, the parameter a decreases: a ~ — ), we have the period 
T (1 THz) ~ v^5- 10 -6 . As is seen, the oscillations of radiation from photonic 
crystals are quite observable. 

So the time r p h = ^- that the photon spends in the crystal can be longer than 

the flight time r p = — of a relativistic particle in a crystal. Hence, the emis- 
sion of diffraction-related radiation (quasi-Cherenkov, transition, diffracted 
radiation of an oscillator, surface parametric radiation and others) produced 
by a relativistic particle will continue after the particle has left the crystal 
(see FigfJ]) Under diffraction conditions, the crystal acts as a high-quality 
resonator [T|l2l]. 

It should be noted, of course, that in observation of oscillations, one should 
either register the moment of particle entrance into the crystal or use a short 
bunch of particles with duration much shorter than the oscillation period. 
In the X-ray range, such situation is typical of electron buches, which are 
applied for creating X-ray FELs (DESY). (The bunch duration in such FELs is 
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tens-hundreds of femptoseconds) . In the terahertz range, much longer bunches 
are required, so there are not serious experimental problems in this case. If 
the bunch duration is large in comparison with the duration of the radiation 
pulse or the time of the electron entrance into the crystal is not registered, 
which occurs in a conventional experimental arrangement, then the intensity 
I(t) should be integrated over longer observation time intervals. As a result, 
we, in fact, obtain the expression (JT]) integrated over all frequencies, i.e., an 
ordinary stationary angular distribution of radiation. If the response time of 
the devices detecting td (or the flight time of the particle in a crystal, or the 
bunch duration) is comparable with the oscillation period, then I(t) should 
be integrated over the interval td- In this case oscillations will disappear, but 
we will observe the power-law decrease in the intensity of radiation from the 
crystal. 

In according with the above analysis some experiments are carried to observe 
delay of radiation pulse in a photonic crystal used for VFEL lasing [251428] . 
In these experiments the additional radiation peak (see Figj2]) is observed at 
studies of lasing of VFEL with "grid" photonic crystals in backward wave 
oscillator regime. This peak appears when the electron beam has left the 
resonator. 




-100 100 200 300 400 500 

Time, nanoseconds 

Figure 2. Detected microwave signal (black curve) synchronized with the beam 
current and electron gun voltage 

It should be mentioned here that backward wave oscillator regime implies 
generation in presence of Bragg diffraction, therefore, under some conditions 
the group velocity could appear even to be close to (see equation ( 135]) ). The 
observed delay (Figj2]) corresponds to v gr ~ 10 8 cm/s, i.e. ^ ~ 1CT 2 . 

In travelling wave regime, which corresponds to case of Laue diffraction, such 
long delay can not be obtained (according to (1351) for (3 > the group velocity 
v gr changes insignificantly). Particularly, in our experiments with Cherenkov 
generator without diffraction grating no additional peaks are detected, because 
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the group velocity in this case changes insignificantly due to the same reasons 
as in the Laue case. 

And after all note that diffraction of a pulse of radiation produced by an 
external radiation source in a periodic structure could be accompanied by 
appearance of several transmitted or reflected radiation pulses (pulses of pho- 
tons) (see [29]). 



Conclusion 

The formulas which describe the time evolution of radiation produced by a rel- 
ativistic particle moving in a crystal are derived. It is shown that the conditions 
are realizable under which parametric (quasi-Cherenkov) radiation, transition 
radiation, diffracted radiation of the oscillator, surface quasi-Cherenkov and 
Smith-Purcell radiation last considerably longer than the time r p of the particle 
flight through the crystal. The results of carried out experiments demonstrate 
the presence of additional radiation peak appearing after the electron beam 
has left the photonic crystal. 
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